BDM Dark Matter: CDM with a core profile and a free streaming scale 
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We present a new dark matter model BDM which is an hybrid between hot dark matter HDM 
and cold dark matter CDM, in which the BDM particles behave as HDM above the energy scale E c 
and as CDM below this scale. Evolution of structure formation is similar to that of CDM model but 
BDM predicts a nonvanishing free streaming Xf s scale and a inner galaxy core radius r core , both 
quantities determined in terms of a single parameter E c , which corresponds to the phase transition 
energy scale of the subjacent elementary particle model. For energies above E c or for a scale factor a 
smaller then a c , with a < a c < a eq , the particles are massless and p redshifts as radiation. However, 
once the energy becomes E < E c or a > a c then the BDM particles acquire a large mass through 
a non perturbative mechanism, as baryons do, and p redshifts as matter with the particles having 



a vanishing velocity. Typical energies are E c = 0(10 — 100)eV giving a A/ s oc E t 



-4/3 < 



Mpc and 



Mf 3 oc i5~ < 10 MQ. A \f s 7^ 0, r core 7^ help to resolve some of the shortcomings of CDM such as 
overabundance substructure in CDM halos and numerical fit to rotation curves in dwarf spheroidal 
and LSB galaxies. Finally, our BDM model and the phase transition scale E c can be derived from 
particle physics. 



I. INTRODUCTION 

The understanding of our universe has received a great 
deal of attention in recent times. Cosmological data 
such as large scale structure [1], SNla [2] and CMB [3] 
are consistent with the concordance ACDM model, with 
fl dmo = 0.22, fl Ao = 74, fl bo = 0.04 with h = 0.71. Even 
though the simple concordance model describes well our 
universe the nature of dark energy DE and dark matter 
DM, which account for up to 96% of the energy content, is 
not well understood. A large number of candidates have 
been proposed for DM of which cold dark matter (CDM) 
has been the most popular. CDM model has been suc- 
cessful on large scales in explaining structure formation 
in the early universe as well as abundances of galaxy clus- 
ters [1]. However, CDM predicts steeply cusped density 
profiles and causing a large fraction of haloes to survive 
as substructure inside larger haloes [4, 5]. These charac- 
teristics of CDM haloes, however, seem to disagree with 
a number of observations. The number of sub-haloes 
around a typical Milky Way galaxy, as identified by satel- 
lite galaxies, is an order of magnitude smaller than pre- 
dicted by CDM [6] and the observed rotation curves for 
dwarf spheriodal dSph and low surface brightness (LSB) 
galaxies seem to indicate that their dark matter haloes 
have constant density cores [7, 8] instead of steep cusps 
as predicted by the NFW profile. Low surface brightness 
galaxies are diffuse, low luminosity systems, with a total 
mass believed to be dominated by their host dark matter 
halos [9]. Assuming that LSB galaxies are in dynami- 
cal equilibrium, the stars act as tracers of the gravita- 
tional potential, and can therefore be used as a probe of 
the dark matter density profile [10]. Much better fits to 
dSph and LSB observations are found when using a cored 
halo model [11]. Cored halos have a mass-density that 
remains at an approximately constant value towards the 
center. 



It has been argued that the sub-structure and core 
problems might be solved in a NFW profile once addi- 
tional baryonic physics are taken into account such as 
reionization and supernova feedback. This feedback may 
help to suppress star formation and to decrease central 
densities in low-mass dark matter haloes [12]. However, 
even tough these processes may help to solve the problem 
with the over-abundance of satellite galaxies, the sugges- 
tion that feedback processes can actually destroy steep 
central cusps seems somewhat contrived in detailed simu- 
lations [13]. Due to these discrepancies in CDM, numer- 
ous alternatives to the CDM paradigm have been pro- 
posed. These include broken scale-invariance [14], warm 
dark matter [15], scalar field dark matter [16], and var- 
ious sorts of self-interacting or annihilating dark matter 
[17]. However, these alternatives are unable to solve both 
problems simultaneously. 

Here we propose a new version of dark matter, well mo- 
tivated from particle physics, which predicts simultane- 
ously a cut in the substructure formation and cored cen- 
ter galaxies. The model simply consist of particles that 
at high energy densities are massless relativistic particles 
with a velocity of light, v = c, but at low densities they 
acquire a large mass, due to nonperturbative quantum 
field effects, and become non relativistic with a vanish- 
ing (small) dispersion velocity. We will name this type 
of dark matter BDM, from bound states dark matter. 
The name is motivated by the particle physics model, 
discussed in section III, but we would like to stress out 
that the cosmological properties of BDM do not depend 
on this particle model but on the different behavior of 
the BDM particles. The phase transition energy density 
is defined p c = and its value can be determined theo- 
retical by the particle physics model or phenomenological 
by consistency with the cosmological data. 

There are two natural places where one encounters 
high energy densities for dark matter. Firstly, at early 
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times when the universe is hot and dense, and secondly in 
the center of galaxies. From a cosmological point of view, 
our BDM has then two clearly distinct behaviors: one 
as hot dark matter HDM at high energies p c < ptdm ~ 
a(i)~ 4 with flbdm being constant and another as CDM for 
energies smaller than p c > pbdm ~ a(<)~ 3 . The transition 
takes place at a c and for a < a c the BDM particles are 
relativistic but for a > a c they become non-relativistic. 
If BDM accounts for all DM than the redshift a c , at 
which we have the phase transition p(a c ) = p ci must 
be smaller than the matter-radiation equality a eq and 
Pc > Peq- From constraints on extra relativistic energy 
density at the time of nucleosynthesis we will have an 
upper limit for a c or equivalently a lower limit for p c . 

Since the particles of BDM travel at the speed of light 
at energy densities above p c there will be a cut at small 
scales in the power spectrum and BDM will erase inhomo- 
geneities and inhibit structure formation for scales below 
the free streaming scale A/ s . This property is similar to 
that of WDM and it predicts lower number of substruc- 
ture in DM halos as a CDM model, welcome by the data. 
The value of A / s will depend only on the phase transition 
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energy scale E c = p c . Once the universe expands and 
Pbdm drops to values smaller than p c , the speed of the 
BDM particles vanishes and these particles will then be 
CDM until present time. 

Once DM dominates the universe, structure formation 
is effective, our BDM particles are cold and we expect 
a standard CDM inhomogeneities growth but with the 
a cutoff A/ s in the power spectrum. The average en- 
ergy densities in halos is of the order 10 5 p o and as long 
as p g < p c we expect a standard CDM galaxy profile, 
which may be given by the NFW profile. However, once 
we approach the center of the galaxy the energy den- 
sity increases and once it reaches the point p g — p c we 
encounter a phase transition for the BDM particles and 
they become massless again. We identify this energy den- 
sity of the galaxy with the core energy density p core = p c 
at a core galaxy radius r = r core . Inside r < r core the 
BDM particles are relativistic and the DM energy density 
inside the radius r core remains constant avoiding a galac- 
tic cusp. Of course we would expect a smooth transition 
region between these two distinct behaviors but the ef- 
fect of considering the thickness of this transition region 
is a small and we will not consider it here. As for the free 
streaming scale A/ s , the size of r core depends only on E c . 

Therefore, our BDM model predicts a free streaming 
scale A/ s , a cut off in the power spectrum and a galaxy 
core radius r core all derived in terms of a single parameter 
E c . The value of E c is determined by the particle physics 
model and constraint by the dark matter properties de- 
rived from cosmology. Furthermore, since E c gives the 
phase transition scale of a particle physics model further 
phenomenological consequences could arise. We present 
in section II the cosmological relevant properties of BDM 
model and in section III the particle physics model giving 
BDM. Finally, we present in section IV our conclusions. 



II. BDM COSMOLOGY 

We consider a standard FRW model universe with dark 
energy and dark matter. The dark matter proposed here 
is motivated in section III from particle physics and it has 
a phase transition at an energy scale E c or equivalently 
at an energy density p c = E^. For energy densities larger 
than p c the particles have a peculiar velocity v = c = 1 
and behaves as relativistic dark matter while for pbdm < 
p c the velocity vanishes and pbdm behaves as cold dark 
matter. We encounter naturally two regions in cosmology 
with a high pbdm ■ The first one is at early times when the 
universe is hot and as long as p&d m (a) > p c (a c ), where a c 
is the transition redshfit, and the second region is in the 
center of galaxies for a radius r < r core with pbdm{r) > 

PiXcore) Pc 

The cosmological evolution of the BDM energy density 
and the peculiar velocity in terms of the scale factor a(t) 
is then given by 

fa, \ 4 

pbdm(a) = pc J > Pc , v = l, a<a c (1) 
Pbdm(a) = pbdmo ^ — J < Pc , v = 0, a>a c (2) 
where the subscript o is at present time and 

p c = E c = pbdm(a c ) (3) 
is the critical or phase transition energy density. In 
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FIG. 1: We show the behaviour of BDM and radiation (solid 
(blue) and dotted (red), respectively) and the transition scale 
a c < a eq where BDM changes its evolution from p oc 1/a 3 to 
p oc 1/a 4 . 

terms of a convenient phase transition scale E c = O(10 — 
100) eV, derived from the particle physics model [18] (c.f. 
sec. Ill), we have 

^^) 4io, >°= 4 * ioi °G^) 4m °/^ 

(4) 

If we require pbdm to account for all dark matter the scale 
factor a c must be smaller than a eq , the matter-radiation 
equivalence scale factor. Of course it also possible that 
Pbdm accounts only partially for DM in which case a c 
could be larger than a eq . We can expressed a c in terms 
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FIG. 2: We show the dependence of Q-bdm as a function of E c . 
tibdm is smaller for larger E c and it is constant for a < a c . 



of the transition energy density as 



Pbdn 



1/3 



3/PQ 



o ^ L bd'mo 
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(5) 



ro is the relativistic en- 
- fib the total matter at 
For a c < a eo , eq.(5) 



"eqi 

0.9 eV" with E eq = 0.5 eV and Qbdm(a) re- 



where we have taken a a = 1, O 
ergy density and f2 mo = Qbdmo 
present time and a e? = f2 ro /fi 
gives £? c > 

mains constant for a < a c since the universe is dominated 
by radiation and pbdm behaves in this region also as a rel- 
ativistic fluid. Therefore, Qbdm(a < a c ) — Qbdm{a c ) = 
Qbdmc is constant and in terms of the transition scale we 
have 

_!_). o.»o5 ( ™) 4/3 
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(6) 



with 3i7 2 = p Tc = p r (a c ) + pbdm{a c ) + pb(a c ), 
^bdmo/^mo = 0.22/0.26 and q - 10" 3 < 1. The amount 
of extra energy density is constraint by nucleosynthesis 
"NS". The upper bound is n extra {MeV) < 0.1 - 0.2 
[20] and from eq.(6) we have then the constraint E c > 
(3.3— 1.8) eV and q = a c /a eg < (0.13 — 0.3), respectively. 
We see that NS constraint allows the transition rcdshift 
a c or E c to be quite close to the matter-radiation equiv- 
alence values but with E c > E eq and a c < a eq . We show 
in fig(2) the evolution of p(a) with the pase transition 
at a c . In fig. (1) we show the dependence of Qbdm as a 
function of E c for a < a c . 

It is well know that a HDM model has a cut in the 
power spectrum and does not allowed structure to form 
below the free streaming scale A/ s containing a mass 
M fs = 4:ir Po (\ fs /2) 3 /3. A large M fs > 10 12 M Q , as 
for neutrinos, is ruled out by structure formation but 
a smaller Mf s ~ 10 9 M© would help to solve the problem 
of having too much substructure in a CDM scenario. The 
free streaming scale A / s is given by 



A/.(t) = 



a(t) I 
Jo 

a(t) 



a(t>) 
" da 



a(t) 



a 2 H H r a r 



(7) 



where we have used dt — da/aH and that the velocity 
of the BDM is given by v = 1 for a < a c and v = for 
a > a c . Using eq.(5) with 3H^(a c ) = p T c(a c ) = pc/^bdmc 
and eq.(7) we get a free streaming scale 
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(8) 
(9) 



and a contained mass within a radius A/ s /2 

^po /XfA 3 fmev 



M fs (E c ) 



A/s 
2 



E r 



5.1 x 10 9 M, 



(10) 
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FIG. 3: We show the dependence of Mf s as a function of E c 
and the larger E c the smaller Mf s . 



0.6 
. 4 
0.2 



-0.2 
-0.4 
-0.6 



20 



40 



80 



100 



(eV) 



FIG. 4: We show the dependence of A/ s as a function of E c 
and the larger E c the smaller A y s . 



We see that A/ s and Mf s only depend on E c and are 
inversely proportional to it. The larger the scale E c (or 
p c ) the smaller Xf s ,Mf s will be. The value E c ~ 60 eV 
gives precisely a M/ s of the order of 10 9 M Q as required 
to inhibit extra substructure in halo but small enough 
to allow for small galaxies to be produced. For scales 
A < Xf s there will be little power on the spectrum. The 
constraint to have Mf s < 10 10 M Q gives a lower bound 
to E c , with E c > 50 eV or p c > (50 eV) 4 . In fig. (3) and 
(4) we show the dependence of Mf s , A/ s as a function of 
E c . 

In a spherical DM distribution the rotation velocity vel 
can be determine straightforward by V 2( f> = Airp with 
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$ = -GM(r)/r, giving 



vel = 



GM(r) 



(11) 



and M — An J p(r) r 2 dr. The Jeans scale radius 
rj(E c ) = \j/2 is given by 



rj(E c ) 



AirGp 



2 7T 

3E 2 



60 eV 

E r . 



11 pc (12) 



where we have used p = p c = E%, the sound speed for a 
relativistic fluid v 2 = 1/3 and 8ttG — 1. The Jeans scale 
implies that inhomogeneitics at scales below Aj inside a 
galaxy are erased by the free streaming of the particles, 
i.e. we have a constant inner galaxy core. The Jeans 
radius is the maximum size for an inner core which co- 
incides with the upper limit for a core radius r < 300 pc 
[9] giving a scale E c > 10 eV. It is well known that in 
a cosmological scenario the Jeans scale r j is larger than 
the Hubble radius dn = 1/H for a < a eq and therefore 
structure formation does not form until a matter domi- 
nation universe. However, in a galactic scenario rj gives 
the largest possible radius for a central core region. In- 
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FIG. 5: We show the inner circular velocity of v as a function 
of r/r core for a LSB type galaxy with p s — 5 x 10 4 p o ,r s — 
3h~ kpc. 
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FIG. 6: We show the circular velocity of v as a function of r 
for a LSB type galaxy with p s = 5x 10 4 p o , r s = 3h~ 1 kpc. 

dependent of the galaxy density profile we can estimate 
the contained cored mass as a function of the core radius 

' core 

using cq.(4) 



M core — 



47rp c r 



3 

core 



(13) 



We see from eq.(13) that a contained mass M core — 
O(1)M requires a small r core = O(l0~ 3 )pc for a E c ~ 
O(10— l00)eV. Since r core is small we expect a small de- 
viation from a CDM profile. There are different galaxy 
profiles in the literature. The standard CDM profile is 
given by the NFW p NFW = p s /[(r/r s )(l + r/r s ) 2 } [4] 
while a core galaxy may be fit with p{r) = p&/[(l + 
r/Vb)(l + (r/rb) 2 )] [7] with pb,rb the core density and 
radius, respectively. The NFW profile has a cuspy inner 
region with pnfw diverging in the center of the galaxy 
with an inner slope a = —1 (p g cx r a ) and p s ,r s are 
galaxy dependent parameters. For typical LSB galaxies 
one has p s ~ 10 4 p o ,r s ~ Zh~ l kpc [4]. As mentioned in 
the introduction, a core inner region with constant p and 
a slope a < —1/2 seems to be preferred by dSph and 
LSB galaxies [9]. The core region is reached at r < r core 
when the galaxy energy density reaches the value p c , 



Pg{^) — Pcoreij'core) — Pc- 



(14) 



The size of the core radius r core depends on the choice of 
the galaxy profile. Since our BDM behaves as CDM for 
p < p c as long as the density of the galaxy is p g < p c we 
expect to have a NFW type profile. Therefore, a possible 
BDM profile is given by a cored CDM profile as 



Pg(r) 



{r CO re/r s + r/r s )(l + r/r s ) 2 



(15) 



with r core <C r s . This profile coincides with pnfw at 
large radius but has a core inner region at r = r core with 
P g (r C ore) = Pc = p s rs/2r C ore giving a core radius 



' core 
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2p c 2£ c 4 ' 
The slope is 

_ dLog[p] r/r cor e(l + 3r/r s + 2r core /r s ) 



a = 



dLog[r 



(1 + r/ ' core )(l + r/r 8 ) 



(16) 



(17) 



and takes the values a — (0, — 1/2, — 1, — 2, — 3) for r = 
(0, r CO re, r cor e < r < r s ,r s ,r s < r) respectively. For 
values of r <C r s , eqs.(ll) and (15) gives in terms of 
r' = rjr C ore a central mass, circular velocity and a slope 

M(r') ~ S7rp c r 3 core - r' + Log[l + r']) (18) 

v(r') * + (19) 

a(r>) * - T ^- 7 . (20) 

Notice that from eqs.(8) and (16) we can express A/ s in 
terms of r core 



A/ s — U I c 0re 



(21) 



with b = (2^ 'p bdmo / 'p s r s ) 1/3 / '(H ^fl ro (l + qj) a propor- 
tionality constant, showing explicitly the interconnection 
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of A/ s and r core . In fig. (5) we show the central circular 
velocity as a function of r/r core while in fig. (6) we plot 
the complete circular velocity as a function of r(kpc). In 
fig. (7) we show the central slope profile as a function of 
r/r core while in fig. (8) we have the complete a as a 
function of r(kpc). 
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FIG. 7: We show the inner slope of the galaxy profile in 
eq.(15) as a function of r/r 
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FIG. 8: We show the complete slope of the galaxy profile in 
eq.(15). 



III. PARTICLE MODEL 

In particle physics there are two different dynamically 
ways to generate a particle mass, namely the Higgs mech- 
anism and the a non-perturbative gauge mechanism. In 
the SM the elementary particles (quarks, electrons, neu- 
trinos) get their mass by the interaction with the Higgs 
field. The dynamically evolution of the Higgs field im- 
plies that at high energies all SM masses vanish but once 
the Higgs settles into the minimum of its potential, at 
the electroweak scale E ew = O(l00GeV), it acquires a 
non vanishing vacuum value giving a mass to the SM 
particles. Therefore, the mass of the fundamental parti- 
cles vanishes at high energies and are non zero below the 
phase transition scale E ew . On the other hand, the non- 
perturbative gauge mechanism is based on the strength of 
gauge interaction. The strength of the coupling constant 
"g" depends on the energy as 



g- 2 (E) = .gr 2 + 8tt 2 6 Log 



E 



(22) 



with gi — g(Ei) and b counts the number of elementary 
particles "Q" charged under the group. If b > 0, as for 
the strong force QCD with SU(N C = 3), N f = 6 (N c , N f 
are the number of colors and flavors, respectively) one 
has b = (llN c — 2Nf)/3 = 7, the gauge coupling constant 
increases with decreasing energy E and we have a non- 
abelian asymptotic free gauge group. The condensation 
or phase transition scale is defined as the energy when 
the coupling constant becomes strong, g(E) 3> I, and 
from eq.(22) we have 



E„ = &e- 8 " 2 / b9 l 



(23) 



The fact that E c is exponentially suppressed compared to 
Ei allows as to understand why E c can be much smaller 
then the initial Ei which may be the Planck, Inflation or 
Unification scale. The scale E c sets a phase transition 
scale where above E c the elementary particles " Q" (e.g. 
quarks in QCD) are (nearly) massless and below E c the 
strong force binds these elementary fields together form- 
ing neutral bounds states such as mesons and baryons 
in QCD. The order of magnitude of the mass of these 
particles is 



tubs = dE c 



(24) 



with d = 0{1) a proportionality constant and a Compton 
wave length 



Xr. 



m BS 



60 eV 

mBS 



2 x 10 _5 cm 



(25) 



In QCD one has E c ~ 200 MeV with the pion mass 
m v ~ lAOMeV while the baryons mass (protons and 
neutrons) nib — 940MeV, i.e. the proportionality con- 
stant is in the range 0.7 < d < 5, and with bound 
mass much larger than the mass of the quarks (m u ~ 
(1 - 3)MeV,m d ~ (3.5 - 6) MeV). Clearly the mass of 
the bound states is not the sum of its elementary particles 
but is due to the non-perturbative effects of the strong 
force and is well parameterized by E c . The dynamical 
formation of bound states is not completely understood 
since it involves non perturbative physics. However, it 
has been shown in RHIC [19] that at high density above 
the transition scale E c , the QCD quarks do indeed be- 
have as free particles, while at low energies there are no 
free elementary quarks and all quarks form gauge neutral 
bound states. Since the interaction strength increases 
at lower energies, the formation of bound states is ex- 
pected to be larger at the smallest possible particle bound 
state energy Ebs (i.e. Ess = mBs) with momentum 
p 2 = E BS — m 2 BS ~ 0. The energy distribution of bound 
states formation is still under investigation [21] and for 
simplicity we take here p = which gives a vanishing 
particle velocity for the bound states. 

It is precisely the non-perturbative gauge mechanism 
that we have in mind for our bound states dark matter 
BDM. Of course, in our case the gauge group and ele- 
mentary fields are not part of the SM. This "dark" gauge 
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group is assumed to interact with the SM only through 
gravity and is widely predicted by extensions of the SM, 
such as brane or string theories. Furthermore, this dark 
gauge group may also account for dark energy as dark 
mesons [18]. The model used in [18] has a transition scale 
E c = O(10 — 100) eV which is the reference energy used 
here. We see from eq.(25) that for m = dE c = 60 eV we 
have a small Compton wave length A com ~ 10~ 5 cto and 
we therefore do not expect it to play a relevant role in 
late cosmology or galactic scales. 

Even tough we have motivated our BDM in terms of a 
well motivated particle physics model we stress the fact 
that the cosmological implications of BDM do not depend 
on its origin of. The BDM is defined by a DM that at 
Pbdm > Pc behaves as relativistic HDM with a particle 
velocity v = l while for pbdm < Pc it is CDM with v ~ 0. 
Naturally one encounters high p at early cosmological 
times and in the inner regions of galaxies. 

IV. SUMMARY AND CONCLUSIONS 

We have presented a new type of dark matter, BDM, 
which can be derived from particle physics. These new 
particles depend on a dark gauge group, similar as the 
SM elementary particles, but with only gravitational in- 
teraction with the SM. The BDM particles acquire a 



non perturbative mass below the phase transition scale 
E c . The order of magnitude for this phase transition is 
E c = O(10 - 100) eV [18]. Above E c we have relativistic 
elementary particles while below this scale we have non- 
relativistic bound states, formed as neutral bound states 
from the elementary fields. The cosmological effect of the 
phase transition manifest itself in two different regions. 
Firstly, in the cosmological evolution of the energy den- 
sity pbdm and secondly in the inner regions of galaxies. 
From a cosmological evolution we have a pbdm redshift- 
ing as radiation above p c = and as matter below this 
scale. Since the BDM particles are relativistic at high 
energies we have a nonvanishing free streaming scale A f s 
with a contained mass Mf s . At the same time our BDM 
predicts a core galaxy profile with radius r core once the 
energy density of the galaxies p g reaches the transition 
scale p c . The quantities Xf s ,Mf s and r core are given by 
eqs.(8), (10) and (16) and have the same origin, namely 
the phase transition scale, and are determined in terms 
of a single parameter E c . Therefore, we have proposed 
a DM model which connects the solution to two of the 
main shortcomings of CDM, namely the overabundance 
of substructure in CDM halos and the rotation curves for 
DM dominated galaxies. 
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